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We propose the magnetic exciton mediated superconductivity occurring in the enigmatic hidden- 
order phase of URu2Si2- The characteristic of the massive collective excitation observed only in the 
hidden-order phase is well reproduced by the antiferro hexadecapole ordering model as the trace of 
the dispersive crystalline-electric-field excitation. The disappearance of the superconductivity in the 
high-pressure antiferro magnetic phase can naturally be understood by the sudden suppression of 
the magnetic-exciton intensity. The analysis of the momentum dependence of the magnetic-exciton 
mode leads to the exotic chiral d-wave singlet pairing in the E g symmetry. The Ising-like magnetic- 
field response of the mode yields the strong anisotropy observed in the upper critical field even for 
the rather isotropic 3-dimensional Fermi surfaces of this compound. 
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Itinerant-localized duality of / electron in condensed 
matter often exhibits an enigmatic phase, whose order 
parameter is almost inaccessible by conventional exper- 
imental techniques [1]. One famous example is given in 
the heavy-fermion superconductor, URu2Si2, where the 
ordered phase signaled by the clear specific-heat anomaly 
at T = 17.5 K has not been resolved over the past 25 
years @]. The magnetic responses of this body-centered 
tetragonal compound are characterized by strong Ising 
anisotropy, and the type-I antiferro magnetic (AFM) 
phase with Q = (1,0,0) emerges in application of pres- 
sure through the lst-order transition from the hidden- 
order (HO) phase jH. 

Interestingly, the bulk superconductivity below T c = 
1.5 K coexists only with the HO phase [4], although sim- 
ilarity of the Fermi surfaces in both the HO and the AFM 
phases is pointed out by Shubnikov-de Haas measure- 
ment Moreover, the upper critical field, H C 2(T) shows 
strong anisotropy, i.e., i? C 2(0) for H | [001] is four times 
smaller than for H || [100], although the 3-dimensional 
Fermi surfaces are rather isotropic [6| . This anisotropy is 
not simply ascribed to theparamagnetic effect due to the 
Ising anisotropy in x(T) [7[, since the NMR Knight shift 
measurement suggests that the dominant contribution of 
the susceptibility is the orbital part and the quasiparticle 
susceptibility is small and rather isotropic [8,0]. 

In the superconducting phase, the power-law T depen- 
dences of the specific heat, C(T) cx T 2 and the 
NMR relaxation rate, Tf 1 cx T 3 at low tempera- 
tures are consistent with the linear density of states at 
low-energy, suggesting the presence of a line node in the 
gap structure. The recent thermal transport [12], |13[ and 
the heat capacity 14J measurements in rotating mag- 
netic fields find the multi-gap behavior and the additional 
point nodes along the c-axis, from which the chiral d- 
wave pairing, k z (k x +ik y ), in the E g symmetry has been 
proposed. 

Recently, the antiferro hexadecapole (AFH) order of 
the U ion has been put forward for the HO phase based 
on the state-of-the-art numerical computation [TBI ]. The 



hexadecapole scenario by the Ginzburg-Landau analysis 
[lrjj and the localized multipole exchange model [ijj have 
provided a natural and coherent description for numerous 
observations. Especially, it provides a convincing expla- 
nation of the collective magnetic excitation observed in 
neutron scattering [l8|, Gil • The mode at Q is observed 
only in the HO phase, but not in the AFM phase (20j as 
a consequence of the strong Ising anisotropy. 

By these circumstances, we immediately suspect that 
the superconductivity is mediated by the massive collec- 
tive magnetic excitation — magnetic exciton. Its realiza- 
tion has been addressed so far in UPd2Al3 2l|, l22j and 
PrOs4Sbi2 23]. In this Letter, we first reproduce the 
characteristic of the observed magnetic-exciton mode on 
the basis of the AFH ordering model. Then, we discuss 
which pairing symmetry is favored by the irreducible- 
symmetry decomposition of the momentum dependence 
of the mode. We demonstrate that the strong anisotropy 
in H C 2 can be interpreted as the Ising-like magnetic-filed 
response of the magnetic-exciton intensity in a coherent 
manner with the AFH- AFM orders. 

Let us begin with the multipole exchange model, 

H t = ~ £ 7(9) [J<p)* q <r- q + D(p)£ q f- q } 
i 

-f5> a) 



where a and £ represent the multipoles proportional to 
the Ising dipole and the xy(x 2 — y 2 )-type electric hexade- 
capole respectively, being active in the assumed low-lying 
r^-IYlf ) (0-50-170 K) crystalline-electric-field (CEF) 
scheme in the last term of (JTJ) [17]. H is the magnetic 
field along the c-axis. Note that J x and J y are inactive 
within this CEF scheme. We assume the linear pressure 
dependences in J(p) and D(p), and the common momen- 
tum dependence 7(9), which is normalized as 7(Q) = 1. 
Although in the previous study the nearest-neighbor in- 
teraction was used in the mean-field (MF) analysis for 
simplicity fl7l ] , here the overall exchange couplings J and 
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FIG. 1. (color online) The calculated magnetic excitation 
spectra at (a) ambient pressure. The symbols are taken from 
[3, (b) p = 0.8 GPa, (c) H/Ho = 0.25, where H is the 
critical field of the AFH phase, and (d) H/Hq = 0.5. The 
strongest intensity at Q in the AFH phase (a) disappears in 
the AFM phase (b). The dashed line indicates the dispersion 
in (a). The increase of H || [001] continuously decreases the 
intensity at Q as in (c) and (d). The spectra in the high-// 
paramagnetic phase (not shown) is similar to (b) in the AFH 
phase. 



D, and 7(g) have been determined to reproduce the ob- 
served phase diagram and the inelastic neutron scattering 
data. To do so, we have required the rather long-range 
exchange interaction up to 7th neighbor, 

7(g) = 2ji{c 2x + c 2y ) + 8j 2 c x c y c z + 4j 3 c 2x c 2y 
+ lGj4C x c y c z (c 2x + c 2y - 1) + 2j 5 (c ix + c iy ) 

+ 4:j e (c ix c 2y + c 2x c iy ) + 2j 7 c 2zi (2) 

where c na — cos(nirq a ) with a = x,y, z. 

The longitudinal dynamical spin susceptibility of f- 
electron, 

/•OO 

X,{q^) = i dt([Sa q (t),Sa. q (0)})e iu \ (3) 
Jo 

with Sa q = (j q — (<jg) can be calculated by using the 
Holstein-Primakoff method at T = 0, EHf. Then, 
the structure function, S(q,uj) — Imx z (9,u)/i, can be 
compared with the observed magnetic-exciton mode in 
the HO phase. Figure Q^a) shows the result of the fitting 
of the inelastic neutron scattering data The overall 
dispersion and the characteristic of the intensity peaked 
at Q and the incommensurate vector Q* = (1, 0.4, 0) are 
well reproduced. The present description of the mag- 
netic excitation spectra provides an interpretation that 
the magnetic exciton is the trace of the CEF excitation 
which propagates by the exchange couplings. 
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FIG. 2. (color online) The MF phase diagrams, (a) p-T and 
(b) p-H . The temperature and the c-axis magnetic field are 
scaled by those of the MF critical values at ambient pressure. 



In the application of the pressure, which increases 
slightly the overall exchange couplings, the property of 
the spectra is gradually changed, and the strongest in- 
tensity at Q in the AFH phase suddenly disappears upon 
entering the AFM phase. As was discussed in [ljj, 
the suppression of the intensity at Q can naturally be 
understood by the change of the role for the transverse 
component at the lst-order phase transition, and it is a 
direct consequence of the strong Ising anisotropy of the 
magnetic moments. 

With the same model parameters, we discuss the 
magnetic-field dependence along the c-axis. The increase 
of H continuously decreases the intensity relatively faster 
at Q than at Q* as shown in Fig.QJc) and (d). Note that 
the dispersion of the shadow band due to the antiferro 
holding gains a slight intensity at H/Hq = 0.5. With 
further increase of H, the spectra in the high-iJ param- 
agnetic phase (not shown) is similar to those of Fig. [ljb) 
in the AFH phase. The tendency of the field dependence 
is consistent with the experiments 26,27]. The obtained 
MF phase diagrams are summarized in Fig. [2l where H 
and T are scaled by the critical values of the AFH phase 
at ambient pressure, Hq and To. 

Now, let us examine the magnetic-exciton mechanism 
of the superconductivity. We consider the effective "du- 
ality" coupling between the heavily renormalized quasi- 
particles and the magnetic exciton as 



Hint = 



-9 2 X* 



(4) 



where X z(q) = Re X z{q,0), s\ = £ Q Y, k ac k a c k+ qa is 
the z-component of the quasiparticle spin density oper- 
ator, and g is the coupling constant between the quasi- 
particle and the magnetic exciton. As the Cooper pair is 
formed by using the attraction mainly in the vicinity of 
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FIG. 3. (color online) The c-axis field dependence of the 
pairing interaction in each symmetry channel. In the sin- 
glet channel and the triplet channel with d(k) || c (d(k) _L c), 
the positive (negative) region means attractive. The E g (yz, 
zx) channel is always dominant in the whole H region. The 
experimental H% 2 is indicated by the vertical dashed line. 



Q, we decompose Xz(q) into the form, 

Xz {k - k') = X z{Q) + xo + X> 7 / 7 (fc)/?(*0, (5) 

7 

where 7 = s,x,y, z,xy,yz, zx,xyz, and the basis func- 
tions of the irreducible representation are defined as 



E g : 
A 2u 

E u : 



c x c v c z 



S X SyC z 



fs(k) = 

fxy(k) 
fz\J^) — C x CyS z ^ 

f X {k) = S xCyC Zl fy{fc) — 



S x CyS z 



(6a) 
(6b) 
(6c) 
(6d) 
(6e) 
(6f) 



Here, c a — cos(nk a ) and s a — sin(7r/c a ), and the irre- 
ducible representations for the odd-parity pairing repre- 
sent the symmetry of the orbital part. The coefficients, 
X<y, are obtained as 

f 1 fjh. r 1 AU' 

Xl = i T I V /7(fe)fc(fe ~ fe,)_X2(Q) ~ Xo]/7(fe/) - 

(7) 

By these decomposition, the pairing interactions for each 
channel 7 are given by Vy = g 2 Xj an d V s — g 2 (Xz{Q) + 
Xo + Xs)- In the singlet channel and the triplet chan- 
nel with d(k) || z (d(k) _L z), the negative (positive) sign 
means attractive because of the Ising form of the interac- 
tion, ((4]). The maximum attractive interaction gives the 
highest T c . 

The result of the decomposition under if||[001] is 
shown in Fig. [3] In the whole region of H , the maxi- 
mum attractive interaction is obtained in the E g sym- 
metry of {yz, zx) type. The degeneracy of an arbitrary 



linear combination of the two-dimensional order param- 
eter in E g at T = T c will be lifted below T c . The stable 
linear combination depends on microscopic details such 
as the topology of the Fermi surface. From the pure ener- 
getics of the condensation energy for the isotropic Fermi 
surface, the chiral pairing is the most favorable 28], 



f yz (k) + if zx (k) cx 



k z c 



(k x + k y )a 



[k x ky)a 



(8) 



which has the line nodes in the basal plane and the 
point nodes along the c-axis. In real space, the Cooper 
pair is formed between the nearest-neighbor interlayers. 
This pairing symmetry is indeed proposed by the thermal 
transport and the heat capacity measurements 12H14 1 . 

Let us examine the anisotropy in the magnetic field de- 
pendence of the pairing interaction. Owing to the Ising- 
like field response of the magnetic exciton, the pairing 
interaction is almost independent of the field in the basal 
plane, while it is weakened by the c-axis field as shown 
in Fig. [3J It is the origin of the strong anisotropy in H c2 
as shown below. 

According to the Pesch approximation, which can be 
used to determine H c2 (T) due to the orbital contribution 
[29l | , the condition of H c2 is given by 



ln|£ 



OO 1 



In{T, H) 



1 +2vrT> ^ — i=0, 



(9) 

where A = p F |V r 7 (0)| and T c = (2u c e~< /ir)e~ 1/x are the 
dimensionless pairing interaction and T c at H = 0. In 
the case of H \\ c, the field dependence of the interaction 
is taken into account in the second term, while it can be 
neglected for H _L c. ui n — nT(2n + 1) is the fermionic 
Matsubara frequency, and 

I n (T,H) = V^(u n W(m n )\ip 7 (k)\ 2 ), u n = 1^-, 

vf sin 

(10) 

with the Faddeeva function, W(z) = e~ z "erfc(— iz) and 
the magnetic length, A = (2|e|£f) -1 / 2 . The bracket 
represents the angular average over the Fermi surface, 
where 9 is the polar angle from the c-axis. We have used 
f{k) = \/l5/4sin(26>) sin0 corresponding to the E g sym- 
metry, and the isotropic Fermi surface with the velocity 
Up for simplicity. 

Figure |4] shows the calculated H c2 anisotropy with the 
experimental ratio, H^ 2 /H = 0.361 @. The reduction 
of the pairing interaction A by H decreases not only the 
absolute values of H c2 and T c (H = 0), but also enhances 
the H c2 anisotropy as in ([9]) . The inset of Fig. [4] shows 
the enhanced anisotropy due to the reduction of A. With 
the moderate coupling constant, A = 0.05, can reproduce 
the observed anisotropy indicated by the symbols Q. 

Ordinary itinerant spin fluctuation is weakly influ- 
enced by the magnetic field because of its itinerant nature 
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FIG. 4. (color online) The calculated H C 2 anisotropy. For 
H || c, the reduction of the pairing interaction as in Fig. [3] is 
taken into account. With the moderate coupling constant, 
A, can reproduce the observed anisotropy indicated by the 
symbols 6]. The inset shows the A dependence of H C 2 for 
H\\ [001]. 



and relatively high energy scale. On the contrary, the 
magnetic exciton is inherent in the localized character, 
and it is rather sensitive to magnetic fields, reflecting eas- 
ily in the superconducting stability via the pairing inter- 
action. The feedback effect due the superconducting gap 
opening may be weak on the massive magnetic exciton 
spectra [3p| . In reality, URu2Si2 is the low-carrier com- 
pensated metal in the HO phase, and hence the multi-gap 
pairing has been discussed 12Hl4| . The extension of the 
present analysis to the multi-gap pairing would not al- 
ter the pairing symmetry and the qualitative aspect, but 
improve the precise form of the gap structure. Due to a 
lack of the detailed information on the Fermi surfaces, it 
is beyond the present work. 

In summary, we have addressed the magnetic exciton 
mediated superconductivity, which presumably emerges 
in the exotic AF hexadecapole ordered phase of URu2Si2- 
Charactering the observed magnetic-exciton mode as the 
massive collective excitation from the AF hexadecapole 
ground state, we have deduced that the most favorable 
pairing is the chiral d-wave singlet, k z {k x + ik y ), in the 
Eg symmetry with both the line and the point nodes. 
Thereby, the Ising-like anisotropic field response of the 
mode results in the anisotropic reduction of the pairing 
interaction, yielding the strong anisotropy in H c i- The 
present picture based on the AF hexadecapole ordering 
provides a coherent description of the enigmatic prop- 
erty of URu2Si2, and it will shed light on further devel- 
opments. 
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